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Aft-ST  R  A  CX-- 


An  electron  energy  equation  is  formulated  for  weakly  ionized  stagnation 

boundary  layers  (degree  of  ionization  in  the  order  of  0.  1  peijcont  or  less) 

over  highly  cooled  surfaces.  Only  the  cases  of  moderately  negative  wall 

potential  with  respect  to  the  plasma,  where  the  ratio  of  the  electron  to 
.  .  . 


ion  current  densities  is  between  about  10,  are  considered.  The 

equation  is  simplified  by  taking  advantage  of  the  fact  that  the  electron 
energy  equation  is  rather  insensitive  to  the  detailed  variation  of  the 
electron  concentration  profile.  A  solution  of  the  simplified  electron 
energy  equation  is  obtained  analytically  for  the  cases  where  the  electron 
temperatures  are  not  in  equilibrium  with  the  neutral  gas  temperatures. 
The  effects  of  nonequilibrium  electron  temperatures  on  the  electrical 
characteristics  of  the  boundary  layer  are  then  analyzed. 

It  was  found  that  the  equilibration  process  of  the  electron  temperatures 
to  the\.eutral  g4s  temperatures  ks  controlled  by  th£«fratio  of  the  charac¬ 
teristic  electron  energy  condj/qtion  time  to  the  Characteristic  tempera¬ 
ture  equilibration  time.  The  characteristic  residence  Name  of  electrons 
has  been yCon-NNsritionally  lfsed  instead  of  the  conduction  time.  / 


The  ion  current  is'^ound  to  be  dependen^on the  electron  tjzmpiqrature  but 
practically  independent  of  the  potential.  / 
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I.  INTRODUCTION 


The  gases  surrounding  hypersonic  vehicles  at  orbital  velocities  and  in  many 
present-day  high  temperature  experimental  facilities  are  frequently  ionized 
up  to  about  0.  1  percent.  An  understanding  of  the  electrical  as  well  as  ther¬ 
mal  interactions  between  such  weakly  ionized  gases  and  a  solid  boundary  is 
important  in  the  study  of  microwave -plasma  interactions  and  in  connection 
with  the  electro-chemical  diagnostics  in  experimental  facilities. 

The  electrical  characteristics  of  the  aerodynamically  important  stagnation 
region  in  weakly  ionized  gases  have  been  previously  studied  by  Talbot  (Ref.  1) 
and  Chung  (Ref.  2).  Talbot  analyzed  the  regime  where  electrical  sheath 
phenomena  are  pronounced  in  a  range  smaller  than  the  mean  free  path  of  the 
gas,  whereas  Chung  considered  the  case  where  electrical  sheath  effects  are 
important  in  a  range  much  greater  than  a  mean  free  path.  From  both  these 
analyses,  it  is  clear  that  the  electron  temperature  is  an  important  factor 
influencing  the  electrical  interaction.  Although  it  is  well-known  that  the  elec¬ 
tron  temperature  in  the  boundary  layer  may  not  be  in  equilibrium  with  the  ion 
and  neutral  gas  temperature,  no  analysis  of  the  nonequilibrium  electron  tem¬ 
perature  profile  is  available  in  the  literature.. 

In  the  present  study,  the  governing  equation  for  the  nonequilibrium  electron 
temperature  profile  is  formulated  for  the  stagnation  boundary  layer  flow  of  a 
weakly  ionized  gas.  With  some  simplifying  approximations  this  equation  is 
solved  essentially  analytically  for  the  case  where  there  is  no  gas-phase 
recombination  of  electrons  and  ions,  i.  e.  ,  the  gas  phase  recombination  of 
electrons  and  ions  is  frozen.  It  is  shown  that  the  gas -phase  recombination  is 
likely  to  be  frozen  in  many  practical  cases,  and  it  may  be  seen  from  the 
solution  that  a  moderate  rate  of  radiative  recombination  does  not  affect  the 
electron  temperature  profiles. 
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The  analysis  is  restricted  to  cases  where  the  gas  density  is  high  enough  to 
permit  a  boundary  layer  approach  and  the  electrical  sheath  effects  are 
pronounced  over  a  range  greater  than  a  mean  free  path..  The  method  pre¬ 
sented  by  Chung  {Ref.  2)  will  therefore  be  used  to  study  the  electrical  inter¬ 
action  between  the  ionized  boundary  layer  and  the  solid  wall  for  the  case  of 
nonequilibrium  electron  temperature  profiles. 


II.  FORMULATION  OF  THE  PROBLEM 


The  governing  equations  and  boundary  conditions  for  a  highly  cooled, 
stagnation  boundary  layer  flow  of  weakly  and  singly  ionized  gas  will  first  be 
formulated.  The  flow  model  is  illustrated  in  Fig.  1. 

A.  GOVERNING  EQUATIONS 


SHOCK 


When  the  degree  of  ionization  is  much 
below  ~  1  percent  (in  the  order  of 
0.  1  percent  or  less),  neutral  gas  flow  — - 

is  practically  unaffected  by  the  charged  — 
particles  present  since  the  neutral- 
charged  particle  collisions  are  rare 
compared  to  the  neutral -neutral 
particle  collisions.  Thus  the  equa¬ 
tions  governing  the  electrical  behavior 
of  the  weakly  ionized  gases  can  be 
decoupled  from  those  governing  the 
neutral  gas  flow.  The  solutions  of  the 

latter  equations  are  well-known  for  highly  cooled  stagnation  boundary  layers 
and  will  not  be  repeated  here.  In  addition,  since  the  molecular  weights  of 
the  ions  and  neutral  gas  particles  are  of  the  same  order  of  magnitude,  the 
ion  temperature  will  be  assumed  to  be  equal  to  the  neutral  gas  temperature 
and,  therefore,  known  in  the  present  analysis. 


Figure  1.  Flow  Model 


•The  governing  equations  for  the  electrical  characteristics  for  two  dimensional 
and  axisymmetric  stagnation  boundary  layers  with  frozen  gas -phase  electron- 
ion  recombination  are  as  follows: 
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Conservation  of  Electrons 


9C  8  C  /  8C  ' 

pu  — r— ~  +  pv  — jr — “  =  -2-  IpD  -  +  pK  E*C 
8x  8y  9y  y  e  9y  r  e  e, 


(2) 
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The  first  three  equations  were  previously  derived  in  Ref.  2.  The  quantities 
and  ejgE*  in  the  energy  equation  are  energy  source  terms  due  to  elastic 
collisions  between  the  electrons  and  other  particles  and  due  to  the  electric 
field,  respectively.  The  term  c  T  [pD  (9C  / 8y)  +  pK  E*C  ]  represents 

p6  6  6  0  0  0 

the  electron  energy  transfer  due  to  the  drift  of  electrons.  The  term 
X.  (9T  /8y)  represents  the  transfer  of  electron  energy  by  conduction.  This 
equation  will  be  discussed  in  greater  detail  in  subsequent  sections. 


< 
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With  the  solutions  of  the  flow  and  neutral  gas  energy  equations  available, 

these  four  equations  remain  to  determine  the  unknown  functions  CX,  C^,  E*, 

and  T  .  The  equations  are  nonlinear  and  basically  coupled  since  D  and  K 
6  6  6 

are  functions  of  T  ,  and  X  ,  W  ,  and  i  depend  on  C  . 

e  e  el  Je  r  e 

B.  BOUNDARY  CONDITIONS 

The  plasma  at  y  -*•  oo  is  assumed  to  be  neutral  and  to  have  a  given  charged 
particle  concentration.  In  addition,  though  it  is  not  necessary,  the  electron 
temperature  is  assumed  to  be  equal  to  the  neutral  gas  temperature  there. 
The  boundary  conditions  at  y  -*■  oo  are  then: 


C.  =  C. 

l  is 


(6) 


M 
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= _ -  c 
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s 

(8) 

Now,  consider  the  gas -solid  interface.  Here  the  gas  layer  is  composed  of 

ions,  electrons,  and  neutral  gas  particles  in  contact  with  the  solid  surface. 

On  arrival  at  the  interface  from  the  boundary  layer,  the  ions  and  electrons 

recombine  catalytically, with  the  wall  itself  acting  as  the  catalyst.  If  the 

surface  is  highly  cooled  and  catalytic  for  electron-ion  recombination,  the 

normalized  electron  and  ion  concentrations  at  the  interface,  C  /C  and 

A  1  A  ew  es 

C.  /C.  „  are  of  the  order  of  (X  / 6)  or  less  ,  where  X  and  6  are  the 
iw  is  gw  g 

neutral  gas  mean  free  path  and  the  boundary  layer  thickness,  respectively. 
For  all  practical  purposes,  therefore,  the  following  wall  boundary  conditions 


1 


This  will  be  shown  in  Sec. 


VI. 
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can  be  used  for  the  conservation  equations  when  the  surface  is  highly  cooled 
and  catalytic: 

At  y  =  0, 


C. 
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CT“ 

is 
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(9) 


C 
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es 

It  can  also  be  shown  (Sec.  VI)  that  for  such  surfaces  the  following  boundary 
condition  is  applicable  for  the  electron  energy  equation: 

At  y  =  0, 


9T 
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"W 


=  0 


(ID 


Boundary  condition  (11)  implies  that  the  temperature  of  the  electrons  at  the 
interface  will  usually  be  different  from  the  surface  temperature.  Thus,  a 
sort  of  "electon  temperature  jump"  exists  at  the  surface.  A  similar 
phenomenon  occurs  in  rarefied  gas  flow  with  regard  to  the  temperature  of 
the  gas  adjacent  to  the  surface. 

The  temperature  jump  becomes  negligible  when  (C ew/Ceg)  >  10  ^ ,  approxi¬ 
mately.  Such  will  be  the  case  if  either  the  surface  is  a  poor  catalyst  for  the 
electron-ion  recombination  or  the  surface  temperature  is  much  higher  than 
1000°K  and  has  a  very  low  work  function  for  thermionic  emission  at  the  same 
time . 

Finally,  Poisson's  equation  is  first  order  and  satisfies  only  one  boundary 
condition.  When  the  surface  is  an  electrically  floating  conductor  or  a  dielec¬ 
tric,  the  boundary  condition  ( j  / j j) w  =  1  can  be  applied  to  the  system  rather 
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than  specifying  E*  at  any  point.  The  solution  in  this  case  will  give  E*  =  0  at 
the  boundary  layer  edge  and  an  E*  corresponding  to  the  Langmuir  (floating) 
potential  at  the  surface.  When  a  net  current  is  flowing,  one  may  give  a  value 
of  E*  or  (j  /j.)  ^  1  at  the  surface.  In  this  case  the  solution  gives  non- 

vanishing  E*  at  the  boundary  layer  edge.  In  Ref.  2  it  was  stated  that  E* 
vanishes  at  the  edge  of  the  boundary  layer.  This  statement  must  be  qualified 
because  E*  vanishes  identically  only  when  E*^  is  that  value  corresponding  to 
( Je / Ji) w  =  ^ •  The  value  of  Eg  and  its  effect  on  the  potential  drop  will  be 
subsequently  discussed. 
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III.  ANALYSIS  OF  ELECTRON  ENERGY  EQUATION 


Before  actually  analyzing  the  electron  energy  equation,  we  should  mention 
that,  due  to  the  considerable  uncertainties  involved  in  evaluating  the  basic 
physical  properties  of  a  plasma,  it  is  not  our  intention  to  obtain  detailed 
numerical  solutions  but  rather  to  derive  an  approximate  but  essentially 
analytical  solution,  permitting  the  study  of  the  effects  of  the  major 
parameters. 

As  explained  in  Sec.  IIA,  the  governing  equations  [Eqs.  (1)  through  (4)]  are 
basically  coupled.  Although  it  is  clear  from  the  Couette  flow  analysis  of 
Ref.  2  that  the  solutions  of  the  first  three  equations  depend  critically  on  the 
electron  temperatures,  the  solution  of  the  electron  energy  equation  is  rather 
insensitive  to  the  detailed  variation  of  the  electron  concentration  profile.  The 
latter  fact  will  be  substantiated  in  Sec.  V.  In  order  to  simplify  the  analysis, 
therefore,  the  electron  concentration  profile  is  considered  to  be  given  by  that 
obtained  from  the  ambipolar  solution  of  the  conservation  equations  [Eqs.  (1) 
and  (2)]  all  the  way  to  the  wall.  From  the  analysis  in  Ref.  2,  it  can  be  seen 
that  such  an  approximation  is  quite  acceptable  for  the  present  cases  where 
the  sheath  is  much  thinner  than  the  boundary  layer.  This  approximation 
decouples  the  electron  energy  equation  from  the  conservation  and  Poisson 
equations,  and  the  analysis  of  the  electron  energy  equation  is  simplified 
considerably.  After  the  electron  temperature  profiles  are  obtained,  the 
electrical  characteristics  of  the  boundary  layer  will  be  studied  by  solving  the 
conservation  and  Poisson  equations  using  a  procedure  similar  to  that  in 
Ref.  2,  with  the  electron  temperature  profiles  found  by  the  method 
described  above. 

An  outline  leading  to  the  solution  of  the  electron  energy  equation  is  given  in 
the  remainder  of  this  section.  The  details  of  the  method  of  solution  are  given 
in  Appendix  A. 


The  present  analysis  is  limited  to  the  regime  in  which  10  *  <  (j'  /'j^)  <  10 

and  (T  /T  )  <  10.  The  former  implies  stagnation  surfaces  whose  potential 
®  § 

with  respect  to  the  plasma  is  weakly  or  moderately  negative.  The  latter 


follows  from  the  already  stated  condition  of  T  =  T  . 

1  es  s 


(Because  T  /T  is  in 
s  w 


the  order  of  10  for  the  highly  cooled  walls,  the  maximum  ratio  of  T  /T  will 

eg 

also  be  of  the  order  of  10.  )  In  this  regime,  the  electron  energy  equation  is 
simplified  substantially  by  observing  the  fact  that  the  energy  source  term 
ejeE*  and  the  electron  energy  transfer  term  due  to  the  electron  drift 
c  T  [pD  (3C  / 9y)  +  pK  E*C  ]  are  negligible  compared  to  the  electron 

p6  6  C  6  0  6 

energy  conduction  term. 


The  electron  thermal  conductivity  is  obtained  from  Fay  (Ref.  3).  The 
expression  for  the  energy  source  term  is  derived  from  the  basic  equa¬ 
tion  given  by  Petchek  and  Byron  (Ref.  4)  with  a  proper  consideration  of  the 
fact  that  the  predominant  elastic  collision  in  the  present  problem  is  between 
the  electrons  and  the  neutral  particles  rather  than  between  the  electrons  and 
the  ions  as  is  the  case  in  Ref.  4.. 


The  electron  energy  equation,  Eq.  (4),  is  thus  simplified  and  becomes 


=  -Pr  r 
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after  the  usual  similarity  transformation, where  , 
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and  (')  denotes  the  differentiation  with  respect  to  the  similarity  variable  T|. 

-3 

For  most  gases,  Pr^  is  of  the  order  of  10  .  The  stream  function  f  is  tab¬ 

ulated  in  Ref.  5.  When  the  collision  cross  section  between  the  electrons  and 
the  air  neutral  particles  given  by  Hansen  (Ref.  6)  is  used,  Eq.  (14)  becomes 
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where  pg  is  in  atmospheres,  Tg  is  in  °K,  and  (3  is  in  sec  \  denotes  the 

ratio  of  the  characteristic  convection  time  (represented  by  1/(3)  to  the  char¬ 
acteristic  equilibration  time  between  the  electron  and  the  neutral  gas  temper¬ 
atures  (represented  by  the  remainder  of  the  terms  in  this  parameter).  On  the 

other  hand,  the  term  Pr  denotes  the  ratio  of  characteristic  conduction  time 

e 

to  characteristic  convection  time  for  the  electron  thermal  energy.  It  is  seen 
that  for  given  h(q)  and  me(t)),  the  solution  of  the  electron  energy  equation 
depends  on  the  two  parameters  Pr  and  T  ..  For  a  given  Pr  ,  T  .  —  0  implies 
that  the  temperature  equilibration  process  is  frozen  compared  to  the  convec¬ 
tion  process,  whereas  r  ^  -*  oo  implies  that  the  electron  temperatures  are  in 
equilibrium  with  the  neutral  gas  temperatures. 

To  solve  Eq.  (12)  analytically,  the  electron  energy  boundary  layer  is  divided 

into  three  regions  as  shown  in  Fig.  2.  The  h  and  mg  profiles  are  each 

replaced  by  linear  profiles.  The  quantity  pg/p  appearing  on  the  right-hand 

side  of  Eq.  (12)  follows  m  instead  of  h  in  region  2  and  is  Pw/pg  *n  region  1. 

The  values  of  q  and  q  that  divide  the  three  regions  are  given  by  Eq.  (25). 
'mm  ‘m 

The  replacement  of  boundary  layer  profiles  by  linear  profiles  is  a  quite 
common  practice,  and  the  m  and  h  profiles  differ  very  little  when  Sc  is  of 
the  order  of  Pr^.  These  approximations,  therefore,  should  not  lead  to  large 
errors  in  the  results.  The  solutions  of  Eq.  (12)  for  the  three  regions 


respectively  are  obtained  analytically  and  satisfy  the  boundary  conditions 

0'  =  0  and  0  =1.  The  results  are: 
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where 
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and  1^  and  are  the  modified  Bessel  and  Hankel  functions,  respectively,  of 
order  v  for  the  variable  p^.  (See  Ref.  7.  )  The  integrals  of  Bessel  and 
Hankel  functions  appearing  in  Eq.  (16)  are  tabulated  in  Ref.  8. 
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where 


A 


=  2 


(19) 


and  I  and  are  the  modified  Bessel  and  Hankel  functions,  respectively,  of 
order  v  for  the  variable  p^. 

For  region  3  (r|  <  p) 


where 
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el 


N  = — ^  (even  integer) 
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(21) 


and  is  the  Hermite  polynomial  given  by 

HN(p3)  =  (2p3)N  -  (N)  (N  -  l)(2p3)N-2 


,  (N)(N  -  1  )(N  -  2)(N  -  3) 

2! 


(2p3) 


N-4 


(22) 


For  extremely  large  values  of  N,  appropriate  asymptotic  forms  of  the 
Hermite  polynomials  can  be  found  in  Ref.  (9)  and  (10). 
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The  four  constants  of  integration,  -y^  through  are  determined  by  requiring 
that  9  and  0'  be  continuous  between  regions  1  and  2,  and  2  and  3,  respectively. 

The  preceding  solutions  were  obtained  for  the  wall  boundary  condition  6^  =  0, 

which  leads  to  an  electron  temperature  jump  at  the  wall.  As  mentioned 

previously,  this  wall  boundary  condition  must  be  replaced  by  the  "no-jump" 

_  2 

condition  of  8  =  h  for  the  special  cases  of  m  >10  .A  solution  of  the 

w  w  ew 

electron  energy  equation  satisfying  the  boundary  condition  8  =  h^  can  be 

readily  obtained  in  the  same  manner  as  the  preceding  ones.  In  the  present 
study,  however,  we  shall  not  analyze  this  special  case  except  for  -*  0. 

The  solution  of  Eq.  ( 12)  for  T  ^  0  and  a  given  9^  is : 

For  regions  1  and  2 


8  = 


(1-8  )  lnll  + - -  nl 

w  \  m  11 
ew 


'  ew  f 


,  n/tt  m1  exp 

ln(m-1)  +  - 

1  ew/  2'JpI 


(pL) 


+  e 


w 


(1  -  erf  Psm) 


(23) 


For  region  3 


9=1- 


(1  -  9W)U  -  erf  p3) 


2  \/Pr  exp(-p,  )  . 

e  r '  3  m/  ,  /  - 1  )  ,  , ,  c  , 

-  ln(m  I  +  (1  -  erf  p  ) 

i—  .  '  ew'  K3m' 


n/tt 


m' 


3m' 


(24) 


The  linear  profiles  for  h  and  m  are  chosen  to  give  correct  values  of  h  and 
r  e  e  w 

h*  ,  and  m  and  m1  ,  respectively.  The  values  of -n  and  it  ,  therefore, 
w  ew  ew  '  'mm  m 

are  given  by 


r|  =  h  ri 
mm  w 'm 


(25) 
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Figure  3.  Nonequilibrium  Electron 
Temperature  Profiles  for  Pr  =  10"^ 


In  this  section,  we  have  analyzed  the  electron  energy  equation,  Eq.  (4),  and 
obtained  the  electron  temperature  profiles  for  the  various  cases.  In  the 
following  section,  the  conservation  and  the  Poisson  equations,  Eqs.  ( 1 ) - ( 3 ) , 
will  be  analyzed  by  using  the  electron  temperature  distributions  that  have 
just  been  obtained. 


IV.  ELECTRICAL  CHARACTERISTICS  OF  BOUNDARY  LAYER 


In  the  present  section,  the  conservation  and  the  Poisson  equations,  Eqs.  (1), 
(2)  and  (3),  will  be  solved.  The  boundary  conditions  to  be  satisfied  are  those 
explained  earlier:  m(oo)  =  1,  m(0)  =  0,  and  the  given  The  electrical 

characteristics  of  the  boundary  layer  will  then  be  studied. 

The  present  study,  like  that  in  Ref.  2,  is  concerned  with  the  regime  in  which 
the  sheath  thickness  is  much  greater  than  the  electron-neutral  particle  mean 
free  path,  but  is  much  smaller  than  the  thickness  of  the  boundary  layer.  The 
method  of  solution  is  essentially  the  same  as  that  given  in  Ref.  2  except  that 
the  more  general  cases  of  unequal  ion  and  electron  temperatures  will  be 
analyzed  here  where  only  the  case  of  equal  ion  and  electron  temperatures  was 
studied  in  Ref.  2. 


The  boundary  layer  is  first  divided  into  two  arbitrary  regions,  the  sheath  and 
the  ambipolar  regions  (see  Fig.  1).  Solutions  will  be  obtained  separately  for 

the  two  regions  and  will  then  be  matched  at  a  suitable  point. 

/ 

In  the  remainder  of  this  section  an  outline  of  the  method  of  solution  and  re¬ 
sults  are  given.  Details  of  analysis  will  be  found  in  Appendix  B. 


Within  the  thin  sheath  region  presently  of  interest,  the  convection  terms  of 

Eqs.  (1)  and  (2)  can  be  neglected.  Also,  the  temperature  ratio 

T./T  =  (T  /T  )  can  be  considered  to  be  constant  in  the  thin  sheath  without 
i  e  g  e' 

causing  any  undue  error.  The  three  coupled  equations,  Eqs.  (1)  through  (3), 
with  the  above  simplifications,  are  integrated  by  using  a  digital  computer  for 
various  combinations  of  controlling  parameters. 


The  ambipolar  region  is  defined  as  that  region  of  the  boundary  layer  where 
the  charge  separation  is  negligible.  Thus,  in  this  region  Poisson's 
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equation  is  extraneous.  The  conservation  equations  (1)  and  (2)  then  give  the 
equation 


(26) 


where 


1 

Sc 

a 


(27) 


and  the  subscript  0  here  represents  the  average  value  for  the  sheath  region. 
For  the  regime  presently  of  interest  in  which  the  ratios  Tg/T^  and  (j  /j1-^) 
are  each  no  greater  than  ~  10,  Eq.  (27)  can  be  simplified  to 


1 

Sc 


(28) 


For  the  case  of  equal  electron  and  ion  temperatures,  the  above  equation  gives 
the  following  familiar  definition: 


Sc  ^  A  Sc.  (29) 

a  2  i 

Consistent  with  the  approximations  that  we  have  made  so  far,  we  consider 
that  Sca  is  constant  across  the  ambipolar  region.  Such  an  approximation  will 
not  cause  any  undue  error  in  the  results  because  the  solution  of  a  diffusion 
equation  such  as  Eq.  (26)  depends  only  weakly  On  the  Schmidt  number 
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(%  Sc^^).  Equation  (26)  then  becomes  the  conventional  boundary  layer 
diffusion  equation  and  the  solution  can  be  readily  obtained  in  terms  of  the 
stream  function  f  (see  Ref.  11). 

For  the  electric  field  intensity  in  terms  of  the  m  profile  obtained  above,  the 
analysis  of  Eqs.  (1)  and  (2)  also  yields  the  following  equation: 


(30) 


where  (T^/T^^  represents  the  average  temperature  ratio  for  the  sheath 
region. 

The  complete  solution  of  the  original  equations,  Eqs.  (1)  through  (3),  for  the 
boundary  layer  is  now  obtained  by  matching  the  numerical  solution  of  the 
sheath  region  and  the  closed-form  solution  of  the  ambipolar  region  at  for 
rm,  m^,  ml,  m^  and  E,  respectively. 

The  potential  drop  across  the  boundary  layer  is  obtained  by  integrating  the 
electric  field  intensity  over  the  physical  thickness  of  the  boundary  layer. 

Figure  4  shows  a  typical  solution  of  the  sheath  region  obtained  by  the  use  of  a 
digital  computer.  Figures  5a  and  5b  show  typical  me>  rm  and  E  profiles  for 
the  boundary  layer.  Figures  6  through  8  show  the  overall  electrical  charac¬ 
teristics  of  stagnation  boundary  layers,  which  will  be  discussed  later. 
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Figure  5b.  Typical  Electron,  Ion,  and  Electric 
Field  Distribution  in  Sheath  Region 


V.  DISCUSSION 


Figure  3  shows  the  electron  temperature  profiles  for  the  various  values  of 
the  parameter  T  ^  where  r^,  defined  by  Eq.  (14),  represents  the  "ratio  of  the 
characteristic  convection  time  to  the  characteristic  temperature  equilibration 
time  as  explained  previously.  As  r  j  is  increased  from  zero,  the  electron 
temperatures  near  the  wall  begin  to  decrease  toward  the  neutral  gas  tem¬ 
peratures  (Fig.  3).  This  initial  relaxation,  however,  causes  the  electron 
temperatures  in  the  outer  portion  of  the  boundary  layer  to  deviate  from  the 
equilibrium  temperatures  (neutral  gas  temperatures).  The  reason  for  this 
is  as  follows.  The  rate  of  temperature  equilibration  for  a  given  small  value 
of  r  ^  is  much  greater  near  the  wall  than  in  the  outer  portion  of  the  boundary 
layer  because  the  gas  density  and  the  differences  between  the  electron  and  the 
neutral  gas  temperatures  are  much  greater  in  the  fluid  near  the  wall  than 
distant  from  it.  The  initial  increase  of  r  j,  therefore,  causes  mainly  tem¬ 
perature  relaxation  near  the  wall.  The  decreased  electron  temperature  near 
the  wall  then  causes  the  temperature  of  the  electrons  in  the  outer  portion  to 
decrease  by  conducting  the  energy  away  from  them. 

This  phenomenon  thickens  the  electron  thermal  boundary  layer  substantially 
(Fig.  3).  The  rather  large  thickening  is  due  to  the  extremely  small  magni¬ 
tude  of  the  electron  Prandtl  number  that  propagates  the  electron  temperature 
disturbances  near  the  wall  through  a  considerable  distance.  Figure  3  shows 
that  the  electron  energy  boundary  layer  remains  much  thicker  than  the  neutral 
gas  thermal  boundary  layer  until  r  ^  becomes  sufficiently  large  so  that  the 
temperature  equilibration  is  nearly  complete  throughout  the  boundary  layer. 
The  complete  temperature  equilibration  is  seen  to  be  accomplished  at  the 
value  of  r  ^  of  ~  10  .  According  to  Eq.  (15),  this  means  that  for  laboratory 
conditions  of  pTg  of  the  order  of  10^  °K /  sec  (commonly  encountered),  the 
complete  temperature  equilibration  is  not  accomplished  unless  the  stagnation 
point  pressure  is  on  the  order  of  100  atm.  These  criteria  for  the  equilibrated 


electron  temperatures  are  seen  to  be  much  more  severe  than  those  commonly 
predicted  previously  (for  instance,  see  Ref.  1).  The  reason  for  this  is  that 
the  conventional  criteria  have  usually  been  based  on  the  comparison  between 
the  characteristic  residence  time  for  the  electrons  and  the  characteristic 
temperature  equilibration  time.  The  characteristic  residence  time  of  elec¬ 
trons  is  related  to  the  electron  current  density  given  by  the  approximate 
expression  of  Eq.  (A-6).  The  minimum  residence  time  occurs  when  (jg/j^)w 
is  maximum  10  in  the  present  cases).  Equation  (A-6)  shows  that  the 
residence  time  presently  of  interest  is  within  about  one  order  of  magnitude  of 
the  ion  diffusion  time.  As  is  mentioned  following  Eqs.  (A-6)  and  (A-7),  the 
parameter  (Pr^/Sc^)  represents  the  ratio  of  the  characteristic  conduction 

time  of  electron  energy  to  the  characteristic  diffusion  time  of  ions;  this 

_3 

parameter  is  of  the  order  of  10  .  Thus,  the  characteristic  conduction  time 

for  electron  energy  is  at  least  100  times  shorter  than  the  characteristic 
electron  residence  time  for  the  ranges  of  ( j g ) w  considered  in  the  present 
study.  Therefore,  it  is  the  electron  energy  conduction  time  rather  than  the 
residence  time  that  should  be  compared  with  the  equilibration  time  in  deter¬ 
mining  the  nonequilibrium  electron  temperatures.  From  these  arguments, 
therefore,  it  is  seen  that  the  parameter  (Pr^^)  that  represents  the  ratio  of 
the  characteristic  conduction  time  to  the  characteristic  equilibration  time  is 
a  more  physically  meaningful  parameter  for  the  prediction  of  temperature 
equilibration  process  than  r  ^  alone.  This  can  be  also  seen  from  the  right- 
hand  side  of  Eq.  (12).  Figure  3  shows  that  the  temperature  equilibration  is 
nearly  complete  when  Pr  r  ^  is  of  the  order  of  10. 

The  semi-broken  line  in  Fig.  3  shows  the  electron  temperature  for  m  =0.1, 

°  r  ew 

0  =  h  ,  and  r  ,  -*•  0  obtained  from  Eqs.  ( 23 )  and  (24).  As  mentioned  earlier, 

w  w  el 

this  is  a  special  case  of  either  a  poorly  catalytic  surface  or  a  highly  seeded 
gas  over  a  high  temperature  surface  in  which  ngw  is  within  about  one  order 
of  magnitude  of  nQS-  These  are  highly  special  cases;  however,  the  figure 
shows  the  drastic  effect  the  wall  could  have  on  the  electron  temperatures  of 
the  boundary  layer  should  such  cases  arise.  Due  to  the  extremely  small 
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value  of  the  electron  Prandtl  number,  the  cooling  effect  of  the  wall  will  be  felt 
by  most  of  the  electrons  in  the  shock  layer. 


It  was  mentioned  earlier  that  the  electron  temperatures  were  expected  to  be 
rather  insensitive  to  the  detailed  variation  of  the  electron  concentration  pro¬ 
files.  This  can  be  substantiated  from  the  electron  energy  equation,  Eq.  (12), 
and  the  solutions  of  this  equation.  The  electron  concentration  profile  enters 
into  Eq.  (12)  through  the  terms  (m^/mg)  and  (p/p,.).  The  density  ratio  was 

made  a  function  of  m  during  the  subsequent  solution  of  the  equation  by  letting 
- 1  6 

p/pg  £  m^  .  The  effect  of  varying  the  electron  concentration  profile  on  the 
solution  of  the  energy  equation  can  be  studied  quite  generally  within  the 
framework  of  the  approximations  given  by  Eqs.  (A-  14)  through  (A- 17)  by 
studying  the  effect  of  varying  mgW  on  the  electron  temperatures.  First,  it 
can  be  seen  from  Eqs.  (A- 14)  and  (A- 15)  that 


0  <  T|  <  T| 

—  —  m 

Tj  >  T) 

'  'm 


as  m  -*■()■ 
ew 


Thus,  m'  /m  is  essentially  independent  of  m'  when  m  -+  0,  which  is  the 
e  e  ew  ew 

present  case.  Next,  the  effect  of  varying  mgW  on  the  right-hand  side  of 
Eq.  (12),  through  its  effect  on  p/pg>  is  approximately  of  the  same  order  as 
the  effect  of  varying  r  It  is  seen  from  the  solutions  of  the  equation  given 
in  Fig.  3  that  a  little  variation  of  r  ^  (such  as  10  to  20  percent)  affects  the 
electron  temperatures  very  little.  Thus,  we  have  now  shown  that  the  electron 
temperatures  are  rather  insensitive  to  the  detailed  variation  of  the  electron 
concentration  profiles.  The  above  arguments  can  also  be  readily  supported 
from  the  analytic  solutions  of  Eq.  (12),  i.e.  ,  Eqs.  ( 1 6) ,  (18),  and  (20),  by 
noting  the  very  weak  dependence  of  the  solutions  on  . 


The  preceding  discussion  leads  to  the  following  observation  concerning  the 
effect  of  the  gas -phase  electron-ion  recombination  on  the  electron 
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temperatures:  According  to  Ref.  12,  it  seems  that  radiative  two-body 
recombination  is  the  predominant  gas-phase  recombination  process  when  it 
occurs  in  gases  such  as  argon.  If  this  is  true,  then  the  effect  of  the  recom¬ 
bination  on  the  electron  temperatures  appears  only  through  its  effect  on  the 
electron  concentration  profile.  It  was  shown  in  the  preceding  paragraph  that 
the  electron  temperature  is  rather  insensitive  to  the  variation  of  electron 
concentration  profile.  Therefore,  a  moderate  rate  of  radiative  gas-phase 
recombination  of  ions  and  electrons  will  not  appreciably  affect  the  electron 
temperature  profile.  If  the  recombination,  however,  is  predominantly  a 
three-body  process  in  which  the  electrons  become  the  third  body  (catalyst), 
another  energy  source  term  must  be  added  to  the  present  electron  energy 
equation  and,  therefore,  the  present  solution  is  invalid. 


Before  leaving  Fig.  3,  we  should  mention  that  the  electron  temperature  pro¬ 
file  for  the  frozen  limit  of  =  0  is  not,  strictly  speaking,  that  of  0  =  1  as 
given  in  the  figure.  In  the  analysis,  the  energy  source  term  ej^E#  appearing 
in  the  energy  equation,  Eq.  (4),  was  neglected  because  it  is  of  much  smaller 
order  of  magnitude  compared  to  the  energy  conduction  term.  This  approxi¬ 


mation,  however,  must  fail  in  the  limit  of  =  0  because  in  this  limit  the 
conduction  term  vanishes  everywhere  with  the  electron  temperature  gradient. 
The  effect  of  the  term  ejeE*  is  to  decrease  the  electron  temperatures  because 
E*  near  the  wall  is  negative,  and  the  electron  current  density  is  toward  the 
wall  for  all  the  cases  considered  in  the  present  analysis.  Therefore,  the 
electron  temperature  profile  in  the  limit  of  -*•  0  will  actually  be  much  like 
that  for  a  small  value  of  T  ^  that  gives  a  small  but  nonvanishing  temperature 
gradient  to  balance  the  ejgE*  term. 


Figures  6,  7,  and  8  show  the  typical  variation  of  the  potential  difference  and 
the  current  densities  with  respect  to  several  major  parameters. 


The  net  potential  difference  across  the  boundary  layer  and  the  current  density 
for  the  electron  temperature  profiles  obtained  in  Fig.  3  and  for  (jg/j^)^  =  l 
assuming  n/M6/M.  =  lO"'  are  shown  in  Fig.  6.  The  potential  difference  is 
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Figure  6.  Floating  Potential  and  Ion 
Current  versus  Equilibration  Parameter 

seen  to  increase  as  the  temperature  equilibration  parameter  F^  is  decreased. 
This  is  to  be  expected,  because  the  electron  temperatures  near  the  wall  in¬ 
creased  with  decreasing  r  ^  in  Fig.  3;  also,  the  sheath  potential  drop  is  seen 
(Ref.  2  and  Table  1)  to  increase  with  the  electron  temperature.  It  is  also 
seen  in  Fig.  6  that  the  ion  current  density  increases  in  a  manner  similar  to 
the  potential  difference  as  the  electron  temperature  increases.  It  is  inter¬ 
esting  to  note  that  a  dependence  of  the  ion  current  density  on  the  electron 
temperature  was  also  suggested  for  the  free-fall  sheath  of  a  static  plasma  in 
the  analysis  of  Ref.  13.  Figure  6  shows  that,  for  pw/pg  between  10  and  20, 
the  floating  potential  of  the  wall  with  respect  to  the  plasma  is  increased  by 
about  4  to  5  fold  as  T  ,  is  varied  from  oo  to  zero.  The  ion  and  electron  cur- 
rent  densities  at  the  same  time  are  seen  to  have  increased  by  about  four  fold. 

Figure  7  shows  the  variations  of  the  floating  potential  and  the  ion  current 
density  with  respect  to  the  plasma  flow  parameter  a^Vg/(31^j  =  |vg/pAgj  for 
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Figure  7.  Floating  Potential  and  Ion 
Current  versus  Plasma  Flow  Parameter 


r  l  =  720  and  =  oo,  respectively.  It  is  seen  that  the  potential  and  the 
nondimensional  ion  current,  as  a  whole,  are  rather  insensitive  to  the  varia¬ 
tion  of  the  flow  parameter  when  the  wall  is  highly  cooled.  These  quantities, 
however,  are  seen  to  vary  considerably  more  with  the  flow  parameter  for 

r  ,  =  720.  The  reason  for  this  is  that,  the  sheath  thickness  varies  with  the 
el 

flow  parameter  and  thus  the  average  temperature  ratio  T^/Tg  for  the  sheath 
is  a  function  of  the  flow  parameter.  This  dependence  of  the  average  temper¬ 
ature  ratio  on  the  sheath  thickness,  when  the  temperature  equilibration  is  not 
complete,  causes  the  potential  difference  and  the  nondimensional  ion  current 
density  to  be  influenced  to  a  greater  extent  by  the  parameter  a-(vg/(31  )  than 
when  T^/T^  =  1.  As  is  seen  in  Fig.  7  and  Ref.  2,  however,  the  potential 

difference  varies  with  the  density  ratio  p  /p  much  more  than  with  the  flow 

w  s 
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Figure  8.  Potential  Across  Boundary  Layer 
and  Ion  Current  versus  Current  Ratio 


parameter.  The  nondimensional  ion  current  density  J3[vg/(P1  )j  is 
independent  of  the  density  ratio  pw/ pg.  The  ion  current  density  will  be 
discussed  subsequently  in  detail. 

Figure  8  shows  the  typical  variation  of  the  potential  difference  across  the 
boundary  layer  and  the  current  density  with  respect  to  the  ratio  (jg/j^-  The 
figure,  however,  maybe  more  conveniently  interpreted  as  showing  the  effect 
on  the  ion  and  electron  current  densities  of  varying  the  potential  difference. 

It  is  seen  that  as  the  potential  difference2  is  increased  from  about  3.  8  to  6.  2 
(for  instance,  at  pw/pg  =  10)  the  nondimensional  ion  current  density  is  in¬ 
creased  only  by  about  10  percent.  The  electron  current  density,  on  the  other 


2Note  that  the  wall  potential  is  negative  with  respect  to  the  plasma. 
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hand,  is  seen  to  have  decreased  by  about  two  orders  of  magnitude,  since  j'^ 
remained  essentially  constant. 

The  following  observations  on  the  current  densities  can  be  made  from  Fig.  6, 
7,  and  8.  The  nondimensional  ion  current  density  g ^ ^ ^ )j  is  a  function 
of  only  the  electron  to  ion  temperature  ratio  and,  thus,  of  the  parameter  F^. 
The  actual  ion  current  density  for  a  given  electron  temperature,  therefore, 
is : 


C. 

jiw^  As»V  M.Sc. 

1  1 

1  /  2 

The  parameter  [P(p  p  )]  is  the  parameter  appearing  in  the  usual  boundary 
s  s 

layer  solutions  for  determining  heat  or  mass  transfer  rate  (see  Refs.  14  and 
15).  C.  is  the  same  as  the  degree  of  ionization  in  the  plasma.  For  a  given 
value  of  electron  temperature  and  [p(p  p  )]  [C.  /(M.Sc.)],  the  variation  of 

S  S  IS  IX 

the  potential  drop  does  not  vary  the  ion  current  appreciably  within  the  present 
orders  of  magnitude.  The  electron  current,  however,  is  quite  drastically 
affected  by  the  potential  drop. 

Equation  (30)  shows  that  the  electric  field  intensity  vanishes  identically  at  the 
edge  of  the  diffusion  boundary  layer  (where  m1  =  0)  only  if'(j  /j,^)  =  1.  For 

all  other  cases,  it  shows  that  a  nonvanishing  field  intensity  exists  ~at  the 
boundary  layer  edge.  This  is  consistent  with  the  boundary  layer  character¬ 
istics  and  also  with  the  physical  problem  itself.  The  supply  of  the  charged 
particles  to  the  boundary  layer  edge  from  the  inviscid  region  is  by  convection 
and  the  drift  due  to  the  electric  field.  The  convection  supplies  both  the  ions 
and  the  electrons  at  the  same  rate.  Thus,  an  electric  field  is  needed  at  the 
boundary  layer  edge  to  satisfy  the  condition  of  the  unequal  current  densities 
when  (jg/j^)w  ^  1-  Equation  (30)  shows  that  Eg  is  positive  when  >  1  ■ 

whereas  it  is  negative  when  (jg/j^)w  <  1. 
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Strictly  speaking,  one  can  consider  in  a  physical  problem  a  solid  surface  in 
contact  with  a  semi-infinite  plasma  only  if  the  surface  is  electrically  floating 
and,  therefore,  ( j g )w  =  !•  The  electric  field  intensity  then  must  exactly 
disappear  at  infinity.  The  steady-state  condition  of  (je/jj.)w  ^  1.  on  the  other 
hand,  can  be  obtained  only  with  the  help  of  an  external  circuit.  This  means 
that  there  must  be  another  electrode  in  contact  with  the  plasma  at  a  finite 
distance  from  the  present  surface.  Therefore,  the  plasma  is  not,  strictly 
speaking,  semi-infinite  and  the  electric  field  intensity  beyond  the  boundary 
layer  edge  is  not  exactly  zero,  although  it  may  be  of  much  smaller  order  of 
magnitude  than  that  within  the  boundary  layer.  Of  course,  there  will  be 
another  boundary  layer  at  the  other  electrode  surface  and  an  appropriate 
potential  difference  across  it. 

Let  us  consider  the  order  of  magnitude  of  the  residual  electric  field  intensity 
and  its  effect  on  the  actual  potential  difference  between  the  wall  and  the 
plasma.  The  magnitude  of  the  electric  field  intensity  at  the  boundary  layer 
edge  relative  to  that  at  the  sheath  edge  can  be  readily  obtained  from  Eq.  (30). 
A  simple  estimate  from  Eq.  (30)  shows  that  the  electric  field  intensity  at  the 
boundary  layer  edge  is  less  than  one  percent  of  that  at  the  sheath  edge  for 
(■’e/ jf)w  -  ^  anc*  ^or  t*le  temPerature  ratios  given  in  Fig.  3.  For  small 
probes  located  at  a  stagnation  point,  the  three  dimensional  effect  will  quickly 
reduce  such  small  electric  field  intensity  to  zero  beyond  the  boundary  layer 
edge;  the  potential  drop  beyond  the  boundary  layer  edge,  compared  to 
that  across  the  boundary  layer,  may  be  neglected.  The  potential  difference 
obtained  for  the  boundary  layer,  therefore,  is  the  potential  difference  between 
the  wall  and  the  plasma,  for  all  practical  purposes,  when  ( j g )w  <  10.  When 
^e^Vw>>10’  h°We  ver,  Eq.  (30)  shows  that  the  residual  electric  field 
intensity  at  the  boundary  layer  edge  can  become  the  same  order  of  magnitude 
as  that  at  the  sheath  edge.  For  these  cases,  the  electric  field  intensities 
that  penetrate  into  the  plasma  in  the  inviscid  region  cannot  be  neglected  and, 
therefore,  the  boundary  layer  treatment  alone  cannot  describe  the  plasma- 
solid  interaction  satisfactorily. 
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VI.  RANGE  IN  WHICH  THE  THEORY  IS  VALID 


In  the  present  section,  the  regime  of  validity  of  the  analysis  based  on  the 
governing  equations,  Eqs.  (I)  through  (4), will  be  examined. 

A.  CONSERVATION  AND  POISSON  EQUATIONS 

There  are  three  criteria  that  must  be  satisfied  in  order  that  the  present 
solution  of  Eqs.  (1)  through  (3)  be  valid: 

1.  Sheath  thickness  must  be  much  greater  than  the  mean  free 
path. 

2.  The  sheath  must  be  sufficiently  thin  compared  to  the 
boundary  layer  so  that  the  convection  effect  can  be 
neglected  within  the  sheath. 

3.  Gas-phase  ion-electron  recombination  must  be  frozen. 

The  above  criteria  are  the  same  as  the  first  three  criteria  given  in  Ref.  2. 
The  fourth  criterion  in  Ref.  2  was  Tg  =  T^,  but  the  present  analysis  removes 
this  restriction.  The  first  two  criteria  are  rather  straight-forward;  the 
following  expression  derived  in  Ref.  2  for  argon,  for  instance,  specifies  the 
criteria: 


where 


s 


and  T  is  in  °K,  I  is  in  cm,  n  is  in  particles/cm  ,  and  p  is  in  atmospheres. 

Nothing  can  be  said  with  any  certainty  concerning  criterion  3  because  the  gas- 
phase  recombination  process  is  not  too  well  understood  for  the  gases  in  the 
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present  flow  regimes.  Assuming  that  the  radiative  gas-phase  recombination 
sets  in  before  the  three  body  processes,  the  frozen  recombination  criterion 
may  be  expressed  as 


c  2 

T,  oyn 

diff  1  s  .  ^  , 

-  ~  -  <  <  1 

1 

chem  Jiw 


(32) 


- 13 

where  the  recombination  coefficient  y  for  argon  is  about  10  to 
- 1 0  3 

10  cm  /sec  (see  Refs.  2  and  12),  and  6  is  the  m-boundary  layer  thickness. 

It  is  seen  from  the  expressions  (31)  and  (32)  that  a  rather  large  regime  of 
flow  conditions  can  be  analyzed  by  the  present  theory.  It  is  shown  in  Ref.  2, 
for  instance,  that  the  three  criteria  are  satisfied  for  the  following  typical 
conditions  obtainable  in  an  argon  arc  tunnel  facility: 


p  =  0.  1  atm 
rs 


i  =  1  cm 


Re  =  100 


T  =  7000°K 
s 

Tq  =  500°K 


1. 4  x  1012  <  n  <<  1.  4  x  1015 
—  es 


It  is  also  estimated  in  Ref.  2  (with  some  misgivings)  that  the  above  plasma 
condition  may  give  equilibrated  electron  and  neutral  gas  temperatures  in  the 
boundary  layer.  The  estimate  was  based  on  the  conventional  comparison  of 


-36- 


the  temperature  equilibration  time  with  the  characteristic  electron  residence 
time.  According  to  the  present  analysis,  however,  the  electron  temperature 
will  be  far  from  being  in  equilibrium  with  the  neutral  gas  temperature  for 
the  above  flow  conditions. 

B.  ELECTRON  ENERGY  EQUATION 

There  are  two  criteria  which  must  be  satisfied  in  order  that  the  continuum 
electron  energy  equation,  Eq.  (4),  be  applicable  for  the  weakly  ionized  gases 
presently  of  interest.  The  first  is  that  the  mean  free  path  for  the  electrons 
between  the  electrons  and  the  neutral  gas  particles  be  much  smaller  than  the 
characteristic  length  of  the  flow.  The  mean  free  path  between  the  electrons 
and  the  neutral  gas  particles  is  of  the  same  order  of  magnitude  as  that 
between  the  neutral  gas  particles.  This  criterion  is,  therefore,  essentially 
the  basic  one  for  the  boundary  layer  flow  itself;  the  conditions  for  the  appli¬ 
cability  of  boundary  layer  theory  have  been  well  established  elsewhere  (for 
instance,  see  Ref.  16). 

In  addition  to  the  above  criterion  for  the  basic  boundary  layer  flow,  the 
electrons  must  have  a  Maxwellian  distribution  such  that  the  electron  tem¬ 
peratures  can  be  defined  across  the  boundary  layer  in  a  continuous  manner. 

Although  the  electron-neutral  particle  collisions  aid  the  attainment  of  a 

3 

Maxwellian  distribution  ,  a  conservative  criterion  for  the  distribution  is 
based  on  the  requirement  that  there  be  a  sufficient  number  of  electron- 
electron  collisions  during  a  suitable  characteristic  time  for  the  electrons  in 
the  boundary  layer.  In  order  to  verify  the  Maxwellian  distributions  of  elec¬ 
trons  between  two  electrodes,  Kerrebrock  compared  (Ref.  18)  the  electron- 
eltectron  collision  time  with  the  electron  residence  time.  As  was  shown  in 
the  preceding  section,  however,  the  characteristic  conduction  time  of 
electron  thermal  energy  is  much  shorter  than  the  characteristic  electron 

_ 

It  is  shown  in  Ref.  17  that  electron-neutral  particle  collisions  alone  in 
an  electric  field  lead  to  an  electron  distribution  which  is  usually  between 
Maxwellian  and  Druynesteyn  distributions. 


residence  time  when  (j  /j.)  <  10.  The  second  criterion  for  the  present 

problem  is  therefore  satisfied  when  the  electron-electron  collision  time  is 
much  shorter  than  the  Characteristic  electron  energy  conduction  time  across 
the  boundary  layer.  The  characteristic  electron-electron  collision  time  xee 
is  obtained  from  Spitz er  (Ref.  19)  as 


(33) 


where  r  e  is  in  sec,  Tg  is  in  °K,  and  n^  is  in  particles/cm^.  The  logarith¬ 
mic  term  of  Eq.  (33)  is  calculated  and  tabulated  in  Ref.  19.  The  conduction 
time  of  electron  thermal  energy  t  is  given  by 
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Here  6  is  the  thickness  of  the  electron  thermal  boundary  layer.  9'  and  p  are 
some  suitable  average  values  for  the  boundary  layer..  The  criterion  is  then 
satisfied  when 


~  <<  1  (35) 

c 

Although  the  normalized  electron  mass  fraction  mg  decreases  monatomically 
from  the  boundary  layer  edge  to  the  wall  (Fig.  5),  the  number  density  ng 
remains  essentially  constant  almost  all  the  way  up  to  the  wall  because 
m  ~p  /p.  This  has  been  noted  in  Ref.  Z0.  Using  n  «n  in  Eq.  (33), 
we  can  show  that  criterion  (35)  is  satisfied  for  most  of  the  flow  conditions  of 
interest  to  us,  such  as  the  example  given  in  Sec.  VI.  A. 
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problem  is  therefore  satisfied  when  the  electron-electron  collision  time  is 
much  shorter  than  the  characteristic  electron  energy  conduction  time  across 
the  boundary  layer.  The  characteristic  electron-electron  collision  time  t 
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where  r  is  in  sec,  Tg  is  in  °K,  and  ng  is  in  particle s/ cm  .  The  logarith¬ 
mic  term  of  Eq.  (33)  is  calculated  and  tabulated  in  Ref.  19.  The  conduction 
time  of  electron  thermal  energy  rc  is  given  by 
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Here  6  is  the  thickness  of  the  electron  thermal  boundary  ayer.  0'  and  p  are 
some  suitable  average  values  for  the  boundary  layer.  The  criterion  is  then 
satisfied  when 


ee 
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(35) 


Although  the  normalized  electron  mass  fraction  me  decreases  monatomically 
from  the  boundary  layer  edge  to  the  wall  (Fig.  5),  the  number  density  ng 
remains  essentially  constant  almost  all  the  way  up  to  the  wall  because 
m  ~p  /p.  This  has  been  noted  in  Ref..  20.  Using  n  «  n  in  Eq.  (33), 
we  can  show  that  criterion  (35)  is  satisfied  for  most  of  the  flow  conditions  of 
interest  to  us,  such  as  the  example  given  in  Sec.  VI.  A. 
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C.  WALL  BOUNDARY  CONDITIONS 

Boundary  conditions  (9),  (10)  and  (11)  will  now  be  derived. 

Consider  the  interaction  between  the  electrons  at  the  interface  w  and  the 
solid  surface  ww.  The  electron  mass  balance  can  be  written  as 

X  M  f(nv^  -  (nv)  1  =  fpD  +  pK  E*C  \  (36) 

4  e  /  ew  V  /  ewwj  y  e  9y  r  e  e  f  w 

When  the  surface  is  highly  cooled  and  catalytic,  the  electron  emission  from 

the  solid  is  negligible  and,  in  addition,  practically  no  electrons  return  to  the 

interface  after  collidine  with  the  surface.  Therefore,  n  v  -*•  0  and 

eww  eww 

Eq.  (36),  with  the  aid  of  Eq.  (A -6),  gives 


where  and  6  are  the  neutral  gas  mean  free  path  and  the  boundary  layer 
thickness,  respectively.  Equation  (37)  shows  that 


A  similar  analysis  shows  that 


(39) 


It  is  seen  from  Eqs.  (38)and  (39)  that  boundary  conditions  (9)  and  (10)  are 
justified  for  all  practical  purposes  when  the  surface  is  highly  cooled  and 
catalytic. 
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An  electron  energy  balance  can  be  written  between  the  interface  and  the 
surface  as 


ill  c 
4  e  pe 


(nvT )  -  (iivt)  =  \  — -  +  c  T  ( pD  — -  +  pK  E*C  jl 
\  /ew  \  /eww  e  3y  pe  e\r  e  Qy  e  el 
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With  the  aidtoEqs.  (36)  and  (A^3),  the  above  equation  becomes 

,(l/2)+er 


iu  (nO)  (T  -  T  )  = 

4  e'  'eww  ew  eww'  Pr  IT  I 
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C  m 
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(40) 


The  above  equation  shows  that  (9Tg/9y)w  ®  for  highly  cooled  and  catalytic 

surfaces  since  n  v  -*■  0.  The  boundary  condition  (11)  is  thus  verified. 
Equation  (40)  shows  that  the  electron  temperature  jump  is  maximum  at  the 
wall  when  the  surface  is  highly  cooled  and  catalytic,  and  thus  no  electrons 
reach  the  interface  from  the  wall. 

Let  us  now  consider  the  other  extreme,  where  no  electron  temperature  jump 

exists  at  the  wall.  As  mentioned  previously,  a  no-jump  condition  exists,  for 

instance,  when  either  the  surface  is  a  poor  catalyst,  the  surface  temperature 

is  much  higher  than  about  1000°K  and  at  the  same  time  the  surface  has  a  very 

low  work  function,  or  when  the  gas  adjacent  to  the  high  temperature  surface 

is  highly  seeded.  For  such  cases,  the  rate  of  electron  transfer  from  the 

interface  to  the  surface  and  that  from  the  surface  to  the  interface  are  of  the 

same  order  of  magnitude,  n  v  ~  n  v  .  The  effect  of  the  wall 
°  eww  eww  ew  ew  g 

boundary  condition  on  the  electron  temperature  profile  is  maximum  when 

r  .  -*0.  Thus,  after  substituting  n  v  for  n  v  and  obtaining 

(9Te/9y)w  from  the  frozen  solution  (23),  Eq.  (40)  can  be  readily  solved  for 

the  order  of  magnitude  of  n  that  would  make  (T  -  T  )  -*■  0.  Such 

ew  ew  eww 

3 

The  electron  temperature  profile  is  almost  independent  of  the  wall  boundary 

condition  when  T  ,  >  700. 

el 
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analysis  shows  that,  for  many  of  the  present  flow  regimes,  rnew 
become  the  order  of  10  ^  or  so  in  order  that  the  electron  temperature  jump 
can  be  neglected.  This  means  that  in  order  to  neglect  the  temperature  jump, 
ngw  must  become  approximately  within  one  order  of  magnitude  of  neg- 


VII.  CONCLUDING  REMARKS 


An  electron  energy  equation  was  formulated  for  weakly  ionized  stagnation 
boundary  layers  over  highly  cooled  and  catalytic  surfaces.  A  solution  was 
obtained  analytically  for  electron  temperature  profiles  that  are  not  in  equi¬ 
librium  with  the  neutral  gas  temperatures.  The  effect  of  the  nonequilibrium 
electron  temperatures  on  the  electrical  characteristics  of  the  boundary  layer 
was  then  analyzed  by  solving  the  conservation  equations  of  the  charged  par¬ 
ticles  and  the  Poisson  equation. 

It  is  found  that, when  the  electron  temperature  is  in  equilibrium  with  the 
neutral  gas  temperature  in  the  inviscid  region,  the  potential  difference  across 
the  boundary  layer  and  the  current  densities  increase  about  four  to  five  times 
as  the  electron  temperatures  in  the  boundary  layer  vary  from  the  completely 
equilibrium  state  to  completely  frozen  state  with  respect  to  the  neutral  gas 
temperatures. 

Though  the  variation  of  the  electron  temperatures  affected  the  ion  current  and 
the  potential  difference  alike,  a  variation  of  the  potential  difference  alone  for 
a  given  electron  temperature  was  found  to  have  a  negligible  effect  on  the  ion 
current.  On  the  other  hand,  the  electron  current  was  found  to  be  drastically 
affected  by  the  variation  of  the  potential  difference. 

The  state  of  the  electron  temperatures  with  respect  to  the  neutral  gas  tem¬ 
peratures  in  the  boundary  layer  was  found  to  be  controlled  by  the  parameter 
(Prerej)  which  represents  the  ratio  of  the  characteristic  electron  energy 
conduction  time  to  the  characteristic  temperature  equilibration  time.  This 
parameter  in  the  regime  presently  of  interest,  |l0 <  (je/ j^)w  <  1  o] ,  was 
found  to  be  about  two  to  three  orders  of  magnitude  smaller  than  the  ratio  of 
the  characteristic  electron  residence  time  to  the  temperature  equilibration 
time.  The  latter  ratio  has  been  conventionally  used  to  determine  the  state  of 
the  electron  temperatures.  The  present  analysis  shows,  therefore,  that  the 
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electron-neutral  gas  temperature  equilibration  is  accomplished  at  much 
higher  pressures,  for  instance,  than  those  previously  expected. 

The  special  cases  of  poorly  catalytic  surfaces  and  the  highly  seeded  gas 
adjacent  to  high  temperature  surfaces  were  considered.  It  was  shown  that 
the  cooling  effect  of  the  wall  can  have  a  drastic  effect  on  the  electron  tem¬ 
peratures  of  the  boundary  layer  should  these  conditions  cause  the  electron 
concentration  at  the  wall  to  become  within  one  order  of  magnitude  of  that  at 
the  inviscid  region. 
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NOMENCLATURE 


A 

a 

b 

C 

c 

pe 

D 

E 

E* 


f 

G1’G2’G3 

H 


I 

1 

3 


j 


K 


K 

1 

k 

L 


parameter  defined  by  Eq.  (B-4) 
parameter  defined  by  Eq.  (B-9) 
constant  defined  by  Eq.  (A- 10) 
mass  fraction 

specific  heat  of  electrons  per  unit  mass  of  electrons 
binary  diffusion  coefficient 

nondimensional  electric  field  intensity  defined  by  Eq.  (B-9) 

electric  field  intensity 

absolute  charge  of  an  electron 

Blasius  stream  function 

parameter  defined  by  Eq.  (B-20) 

nondimensional  electric  field  intensity  defined  by  Eq.  (B-4) 

Hermite  polynomial  of  order  N 

T  ,/T 
g  s 

modified  Bessel  function  of  order  v 

nondimensional  ion  current  density  defined  by  Eq.  (B-4) 
nondimensional  ion  current  density  defined  by  Eq.  (B-12) 
particle  current  density 
mobility 

modified  Hankel  function  of  order  v 

Boltzmann  constant 

physical  thickness  of  sheath  region 


-47- 


t 

M 

m 

N 


P 

Pr 

e 

Pr 

Q 


g 


Q 


eg 


Re 

r 

s 

Sc 

a 

Sc 

€ 

Sc. 

1 

T 

u 

V 

V 

A 

V_  . 

je_ 


x 


nose  radius 
particle  mass 
C/C 

s 

even  integer  defined  by  Eq.  (21) 
charged  particle  number  density 
pressure 

electron  Prandtl  number  defined  by  Eq.  (13) 

Prandtl  number  of  neutral  gas 
Coulomb  collision  cross  section 

collision  cross  section  between  electrons  and  neutral 
particles 

Reynold's  number 

distance  shown  in  Fig.  1 

variable  defined  by  Eq.  (A-9) 

ambipolar  Schmidt  number  defined  by  Eq.  (27) 

electron  Schmidt  number 

ion  Schmidt  number 

temperature 

x-component  of  velocity 

potential 

y-component  of  velocity 

electron  mean  thermal  speed  . . 

energy  transfer  to  electrons  per  unit  time  per  unit  volume 
due  to  elastic  collision 

streamwise  distance 
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o 


distance  normal  to  surface 


variable  defined  by  Eq.  (B-4) 


<1-  y2  -  • 


du  /  dx 
s 

temperature  equilibration  parameter  defined  by  Eq.  (14) 
constants  of  integration 

constant  portion  of  Q  defined  by  Eq.  (A-5) 

eg 

modified  thickness  of  sheath  region  defined  by  Eq.  (B-4) 


0  for  two  dimensional  body 
1  for  axisymmetric  body 


variable  defined  by  Eq.  (A- 10) 

T  /  T 

e  s  _ 

Hebye  shielding  length,  "\J (kT^)/(e2ne) 

thermal  conductivity  of  electrons 

viscosity  and  kinematic  viscosity,  respectively. 

density 


Pl’ P2’ P2’ P3 


variables  defined  by  Eqs.  (17),  (A-23),  (19)  and  (21), 
respectively 


SUBSCRIPTS 


characteristic  time 


electrons 


neutral  gas 


matching  point  between  regions  2  and  3 


matching  point  between  regions  1  and  2 


t 


s  edge  of  boundary  layer 

w  wall 

ww  solid  surface 

v  integer 

0  matching  point  between  sheath  and  ambipolar  regions 

unless  specified  otherwise 


NOTE:  Symbols  without  the  subscript  "e"  or  "i"  refer  to  neutral  gas. 


APPENDIX  A 


Supplement  to  Analysis  of  Electron  Energy  Equation 


The  details  leading  to  the  solution  of  the  electron  energy  equation  are  given  in 
the  following  paragraphs. 


I.  Fundamental  Property  Values  for  the  Electron  Energy 
Equation 


The  electron  thermal  conductivity  Xg  and  the  energy  source  terms  W  j  and 
ejfiE*  will  be  studied  in  this  section. 


Consider  first  the  electron  thermal  conductivity  Xg.  Fay  (Ref.  3)'  gives 
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The  Coulomb  cross  section  Qc  is  about  two  orders  of  magnitude  greater  than 

the  electron -neutral  cross  section  Q  so  that  the  second  term  in  the 

e8  _2 

denominator  dominates  for  highly  ionized  plasma  (n  /n  »  10  )  and  X  is 

eg  _2  e 

independent  of  n  .  In  the  weakly  ionized  plasmas  (n  /n  «  10  )  of  current 

e  eg 

interest,  the  first  term  predominates  and  Xg  is  practically  a  linear  function  of 

n  .  It  can  be  seen  from  data  in  Ref.  6  and  other  available  data  that 
e 
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where  -1/2  ^  <r  1/2,  depending  on  the  gas  and  temperature  range.  Equation 
(A-T)  can  be  written  as 


where 


Pr  = 
e 


/  2 


(A-l) 


(13) 


We  shall  call  Pr  ,  defined  by  Eq.  (13),  the  electron  Prandtl  number.  As  will 
e 

be  seen  subsequently,  Prg  influences  the  behavior  of  the  electron  thermal 

boundary  layer  in  a  manner  similar  to  the  usual  influence  of  the  Prandtl 

number  on  the  conventional  thermal  boundary  layer.  The  electron  Prandtl 

-3 

number  is  of  the  order  of  10  for  most  gases. 

As  in  Ref.  2,  the  present  study  is  limited  to  the  regime  10  1  $  ^e^i\v  ^  *0, 
Certain  simplifications  of  the  electron  energy  equation  are  possible  in  this 
regime.  The  electron  energy  source  term  eJeE*  represents  the  rate  at  which 
electrons  gain  and  lose  energy  due  to  their  motion  in  the  electric  field.  Using 
values  of  j  and  E*  obtained  from  Eq.  (4)  and  Ref.  2,  it  can  be  shown  that 
this  term  is  negligible  compared  to  the  conduction  terms  in  the  ambipolar 
region.  The  electron  temperatures ,  therefore,  are  practically  unaffected  by 
the  electric  field  in  most  of  the  boundary  layer.  In  the  sheath,  however,  this 
term  can  become  of  the  order  of,  or  even  greater  than,  the  conduction  term, 
depending  on  the  magnitude  of  E*.  However,  for  10  *  <  (je/jj)w  —  10,  it  can 
be  seen  from  Ref.  2  that  E*  will  be  small  enough  so  that  this  term,  when 
compared  to  the  conduction  term,  can  be  neglected  even  in  the  sheath.  The 
main  reason  for  this  is  that  the  conduction  of  energy  takes  place  very 
efficiently  because  of  the  high  mobility  of  the  electrons  .  This  fact  can  also 
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be  seen  from  the  extremely  small  value  of  Prg  as  given  by  Eq.  (13).  In  this 
analysis,  therefore,  the  term  ejeE*  is  neglected. 

Finally,  the  energy  source  term  is  considered.  This  term. represents 
the  rate  of  energy  exchange  between  electrons  and  the  rest  of  the  gas  due  to. 
elastic  collisions.  Petchek  and  Byron  have  given  (Ref.  4)  a  basic  expression 
for  W  j  although  the  final  expression  they  derived  is  not  applicable  to  the 
present  problem  because  they  were  interested  in  plasmas  in  which  the 
Coulomb  interaction  predominates,  i.e.  ,  in  which  most  of  the  energy 
exchange  is  between  electrons  and  ions.  For  the  weakly  ionized  gases  of 
current  interest,  most  of  the  energy  exchange  is  between  electrons  and 
neutrals.  The  basic  equation  in  Ref.  4,  however,  may  be  used  to  derive  the 
following  equation,  which  is  valid  for  weakly  ionized  plasmas  (much  less  than 
one  percent): 

Wel  *  3^^gV^-(ra^)”eng(Tg  '  V  |A-4’ 

He  e  it  has  been  considered  that  the  collision  cross  section  between  electrons 
and  neutrals  has  the  form: 
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The  above  expression  can  be  derived  from  information  in  Ref.  6  for  high 
temperature  air. 
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II.  Simplification  and  Transformation  of  the  Electron 
Energy  Equation 


The  equation  has  already  been  simplified  by  decoupling  it  from  the  conserva¬ 
tion  and  Poisson  equations  (see  text)  and  by  neglecting  the  term  ejeE*.  One 
additional  basic  simplification  will  be  made  before  attempting  the  solution. 


It  can  be  readily  shown  from  Eqs .  (1)  and  (2)  that  for  the  predominantly 
ambipolar  boundary  layer  considered  here: 


9  C  /j\/  TV  ac 


(A-6) 


This  approximate  equality  is  sufficient  for  the  immediate  purpose  of  showing 
that  the  term  [pD  (9C  /9y)  +  pK  E#C  ]  can  be  neglected  in  the  electron 
energy  equation.  With  the  aid  of  Eqs.  (5),  (A-3),  (13)  and  (A-6)  the  energy 
equation  may  be  written: 
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The  present  analysis  is  concerned  with  cases  in  which  Tfi  and  T^,  and  je  and 
j^,  differ  by  no  more  than  an  order  of  magnitude  or  so  (see  text).  For  these 
cases  the  term  [(j  /j.)  (1  +  T  /T.)(Pr  J/(T  /To)(1/2)+<r(Sc.)]  «  1  and  may 
be  neglected.  This  is  because  the  ratio  Pre/Sc^,  which  may  be  interpreted 
as  an  electron  Lewis  number,  is  much  smaller  than  unity  and,  therefore,  the 
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electron  energy  transfer  (of  the  energy  CpgTe  only)  due  to  concentration 

gradient  and  electric  field  motivated  diffusion  is  negligible  compared  to  the 

energy  transfer  due  to  thermal  conduction.  Such  is  not  the  case,  for 

instance,  if  the  wall  is  strongly  and  positively  biased  with  respect  to  the 

plasma;  then,  (j  /j.)  »  10.  The  electron  energy  equation,  Eq.  (4), now 

C  1  w 

takes  the  form: 
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To  solve  the  above  equation,  we  will  first  transform  the  equation  to  a 
similarity  form  at  the  stagnation  region.  Let 
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where 


b  = 


^s^s 


=  constant 


Also,  we  define  the  following  nondimens ional  quantities: 

9  =  T  /T  ,  h  =  T  /T  ,  m  -C/C  (A- 1.1) 
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With  the  aid  of  Eq.  (A-4)  and  Eqs.  (A-9)  through  (A-ll),  Eq.  (A-8)  then 
become  s 
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(A- 12) 


where  r  ^  is  given  by  Eq.  (14).  The  stream  function  f,  appearing  in 
Eq.  (A- 12),  is  the  solution  of  the  Blasius  momentum  equation: 

f'"  +  ff"  =  0  (A- 13) 

having  the  boundary  conditions  f(0)  =  f'(0)  =  0  and  f 1  (co)  =  2.  The  solution 
of  Eq.  (A-13)  is  found  in  Ref.  5. 


III.  Solution  of  Electron  Energy  Equation 


To  solve  the  equation  analytically,  we  make  the  following  simplifications. 
First,  the  quantity  m  (Q/h)^^+cr  represents  essentially  the  variation  of  X  . 
For  the  values  of  -1/2  <  or  <  1/2,  (0/h)  ^  varies  very  slowly  compared 

to  m  ,  which  varies  from  zero  to  1  across  the  boundary  layer.  It  is  assumed, 
therefore,  that  (0/h)  °"  =  1  when  compared  to  the  variation  of  me.  With 

this  assumption,  Eq.  (A- 12)  is  linearized  as 


0"  +  (—  +  Pr  fV  =  -Pr  F  ,  W  -  0 ) 

\me  e  /  e  el\PS  / 


(12) 


Next,  consider  the  parenthesized  quantity  in  the  second  term  of  Eq.  (12).  It 
is  seen  that  Pfgf  «  m^/me  for  most  portions  of  the  boundary  layer.  Only 
in  the  outer  portion  of  the  boundary  layer  does  Pref  become  of  the  same 
order  as  m^/me,  and,  in  the  limit  of  y  -*-oo,  Prgf  becomes  greater  than 
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m^/m^  as  -*•  0.  Also,  it  was  stated  earlier  that  the  electron  temperature 
profile  will  be  rather  insensitive  to  the  details  of  variation  of  me  profile. 

With  the  above  facts  in  mind,  we  subdivide  the  electron  thermal  boundary 
layer  into  three  regions  (Fig.  2)  and  the  following  simplifying  approximations 
are  made  for  the  respective  regions: 


For  regions  1  and  2  (0  <  r|  <  r|  ) 
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(A- 14) 


(A-15) 


For  region  1  (0  <  r|  <  t^^) 


P  =  P 
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(A-16) 


For  region  2  (ri  <  ri  <  n  ) 
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Ps/P  =  (A- 17) 

ri  and  Tlrnm  are  given  by  Eq.  (25).  The  above  approximations  are 
comprised  of  replacing  the  h  and  mg  profiles  by  a  pair  of  linear  profiles 
and  of  letting  p  / p  essentially  follow  m  instead  of  h.  A  replacement  of 

S  6 
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boundary  layer  profiles  by  linear  profiles  is  quite  conventional  practice,  and 
the  me  and  h  profiles  differ  very  little.  The  approximations,  therefore, 
given  by  Eqs.  (A-14)  through  (A-17)  should  create  no  great  error  in  the 
results.  The  quantities  h'  and  m'  will  be  obtained  from  the  exact  solution 
of  the  neutral  gas  energy  equation  and  the  exact  ambipolar  solution  of  the 
conservation  equations,  respectively.  The  expression  for  the  stream 
function  f  given  in  Eq.  (A- 15)  is  obtained  from  Ref.  11  for  the  outer  region 
of  the  boundary  layer.  As  was  explained  earlier,  mgw  =  0  for  all  practical 
purposes  when  the  surface  is  highly  cooled  and  catalytic.  However, 
will  be  retained  in  the  analysis  for  the  purpose  of  generality. 

It  is  readily  seen  from  Eq.  (12)  that  in  the  outer  portion  of  the  boundary 
layer,  where  m^/m^  «  Pr  f,  Eq.  (12)  becomes  identical  in  form  to  the 
conventional  thermal  boundary  layer  equations  with  a  Prandtl  number  equal 
to  Prg.  Thus,  Prg  influences  the  electron  energy  boundary  layer  in  a  manner 
similar  to  the  influence  of  the  regular  Prandtl  number  on  the  conventional 
thermal  boundary  layer. 

In  the  following,  the  electron  energy  equation,  Eq.  (12),  will  be  solved 
subject  to  the  simplifications  mentioned  above  for  each  region;  the  solutions 
will  then  be  matched  at  hmm  an^  Tlrn  for  0  and  O',  respectively.  The 
boundary  conditions  for  Eq.  (12)  are 

9'(0)  =  0  (A- 18) 

0(oo)  =  1  (A-19) 


For  region  1,  Eq.  (12)  becomes; 
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Equation  (A-20)  is  then  transformed  to  the  following  inhomogeneous  modified 
Bessel  equation  of  zeroth  order: 
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The  solution  of  the  above  equation  that  satisfies  the  boundary  condition 
(A- 18)  is  obtained  as  Eq.  (16). 

jFor  region  2,  Eq.  (12)  becomes: 


Equation  (A- 23)  is  then  transformed  into 
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The  above  is  another  modified  form  of  an  inhomogeneous  Bessel  equation 
and  the  solution  is  given  by  Eq.  (18). 

For  region  3,  Eq.  ( 12),  becomes  the  following  homogeneous  Herrnite 
equation: 


dp  3  K3 


el 


where 

p3  =  \fPre  (q  -  0.  86038) 


(A-25) 
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If  we  let 


where  N  is  an  even  integer,  a  particular  solution  of  Eq.  (A-25)  is  the 
Herrnite  polynomial  given  by  Eq.  (22).  The  general  solution  of  Eq.  (A-25) 
that  satisfies  boundary  condition  (A- 19)  is  obtained  as  Eq.  (20). 

The  constants  of  integration  through  are  determined  by  requiring 
that  0  and  O'  be  continuous  between  regions  1  and  2,  and  2  and  3, 
respectively. 
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APPENDIX  B 


Supplement  to  Electrical  Characteristics  of  Boundary  Layer 


Details  of  the  solutions  of  Eqs.  (1),  (2),  and  (3)  will  be  given  in  the  following 
paragraphs.  As  explained  in  the  text,  the  solutions  will  be  obtained  for  the 
sheath  and  the  ambipolar  regions  (Fig.  1),  and  the  two  solutions  will  be 
matched  at  a  suitable  point.  The  electron  temperature  distribution  is 
assumed  to  be  known  (Appendix  A). 


I.  Sheath  Region 


The  convection  terms  of  Eqs.  (1)  and  (2)  can  be  neglected  in  the  thin  sheaths 
presently  of  interest.  Equations  (1)  through  (3)  then  become  the  equations 
as  derived  in  Ref.  2. 
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where 


a  =  C/C0  (B-4) 


The  boundary  conditions  are  as  follows 
At  z  =  0 

a.  =  a  ~  0 
i  e 

At  z  =  1 

a.  =  ae  =  1  (B-5) 
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In  the  derivation  of  Eqs.  (B-l),  (B-2)  and  (B-3),  it  is  assumed  that 
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Solutions  of  the  above  boundary  value  problem  were  obtained  numerically 
mainly  for  (T./T  )  =  1  by  Chung  (Ref.  2).  To  avoid  excessive  numerical 
complications  and  to  utilize  the  basic  digital  computer  program  vised  by 
Chung, we  assume  that  the  temperature  ratio  T./Tg  is  constant  and  is  equal 
to  the  average  value  within  the  sheath.  Such  an  approximation  should  be 
acceptable  for  the  thin  sheaths  presently  of  interest.  Solutions  of  Eqs.  (B-l) 
through  (B-3)  are  obtained  for  A  =  3000  and  T^/Tg  equal  to  1,  1/2,  1/3,  1/4, 
and  1/6,  by  using  a  digital  computer.  Major  results  are  given  in  Table  1. 
The  reason  for  choosing  this  particular  value  of  A  will  be  explained 
subsequently. 

Although  considerable  analyses  have  been  made  of  the  high  pressure  sheaths 
since  the  time  of  Langmuir,  most  of  them  lack  generality  and  some  of  them 
lack  even  theoretical  soundness.  A  more  detailed  evaluation  of  some  of 
them  is  given  by  Su  and  Lam  (Ref.  21).  To  study  the  electrical  interaction 
between  the  plasma  considered  here  and  a  solid,  we  must  analyze  equations 
of  the  type  given  here  (Eqs.  (1)  through  (3))  with  the  proper  two-point 
boundary  conditions.  The  equations  are  reduced  to  the  type  given  by  Eqs. 
(B-l)  through  (B-3)  for  static  plasmas.  Only  a  few  of  the  existing  analyses 
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consider  the  general  solution  of  these  equations  for  their  respective 
problems.  They  are:  the  static  electric  discharge  problem  (Allis  and  Rose, 
Ref.  22),  the  plasma-solid  interaction  at  the  shock  tube  end-wall  (Sturtevant, 
Ref.  23),  the  interaction  between  a  static  plasma  and  a  spherical  probe 
(Su  and  Lam,  Ref.  21,  and  Cohen,  Ref.  24),  and  the  Couette  and  stagnation 
boundary  layer  problems  (Chung  Refs.  2  and  25).  Turcotte  (Ref.  26)  also 
obtained  a  solution  of  the  equations  similar  to  Eqs.  (B-l)  through  (B-3)  in 
connection' with  a  shock  tube  experiment.  It  is  also  pointed  out  here  that  all 
these  referenced  analyses  used  the  same  wall  boundary  conditions  as  (B-5). 


II.  Ambipolar  Region 


The  ambipolar  region  is  defined  as  that  region  of  the  boundary  layer  where 
the  charge  separation  is  negligible.  In  this  region,  therefore,  rm  =  =  m, 

The  Poisson  equation  is  then  extraneous  and  the  two  conservation  equations 
describe  the  region  in  the  manner  described  below. 


Subtracting  Eq.  (2)  from  Eq.  (1)  and  integrating  once,  we  obtain  at  the 
stagnation  region: 
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where  is  the  constant  of  integration.  By  using  equations  (A-10)  and 
(B-6),  Eq.  (B-7)  can  be  written  in  the  following  form  for  E: 
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where 
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(B-9) 


Here,  the  integration  constant  is  important  because  it  determines,  for 
instance,  the  residual  electric  field  intensity  at  the  boundary  layer  edge 
where  m'  =  0.  As  mentioned  earlier,  the  one  boundary  condition  due  to  the 
Poisson  equation  is  specified  when  the  current  density  ratio  Ue/.i^)w  is 
specified.  We  shall,  therefore,  determine  the  integration  constant  y,-  in 
terms  of  (j g  / j^) w .  Consistent  with  the  preceding  assumption  of  the  negligible 
convection  effect  in  the  sheath  region,  Eqs.  (1)  and  (2)  as  they  approach  the 
sheath  edge  from  the  ambipolar  region  become: 
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( B  - 1 2.) 


With  a  little  manipulation  of  the  above  equations,  we  have  at  the  sheath  edge: 


aE0m0 


(B  - 13) 


After  some  manipulation,  Eqs.  (B-8)  and  (B-13)  give: 


Upon  eliminating  from  Eqs.  (B-8)  and  (B-14),  the  quantity  (Em)  is 
obtained  from  equation  (B-8)  and  is  substituted  into  Eq.  (1).  After  the 
similarity  transformation  from  (x,  y)  to  q,  this  equation  becomes  Eq.  (26). 
The  solution  of  Eq.  (26)  that  satisfies  the  boundary  conditions  m(qp)  =  m^ 
and  m(oo)  =  1  is  shown  in  Ref.  2  as 
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The  slight  difference  in  constants  between  equation  (B-15)  and  its  correspond¬ 
ing  equation  in  Ref.  2  is  due  to  the  difference  in  the  definition  of  f  and  q  of 
the  two  analyses.  The  quantities  m^  and  q q  are  unknowns  and  will  be 
determined  through  the  subsequent  matching  process. 
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Once  the  m-profile  is  obtained,  the  electric  field  intensities  in  the  ambipolar 
region  are  obtained  form  Eq.  (30),  which  is  derived  by  combining  Eqs.  (B-8) 
and  (B-14) . 

III.  Matching  of  the  Two  Regions 

The  detailed  matching  procedure  is  given  in  Ref.  2.  Here  we  will  only 
briefly  explain  the  basis  upon  which  the  matching  is  accomplished.  In  the 
numerical  solution  of  Eqs.  (B-l)  through  (B-3)  it  was  seen  that  the  region 
in  which  the  charge  separation  and  the  electrical  field  intensity  are 
pronounced  becomes  confined  to  that  portion  of  the  sheath  region  adjacent 
to  the  wall  as  A  is  increased  (Fig.  4  and  also  Ref.  2).  The  charge 
separation  and  the  electric  field  intensity  become  very  small,  in  the  rest 
of  the  region.  The  solutions  for  large  A  thus  produce  a  sort  of  ambipolar 
tail  (Fig.  4).  The  following  equations  can  be  derived  from  Eqs.  (B-l) 
and  (B-2)  for  the  ambipolar  tail: 


(B-1.7) 
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The  first  step  in  matching  the  two  regions  is  to  obtain  a  solution  of  Eqs.  (B-l) 
through  (B-3)  that  has  a  sufficiently  long  ambipolar  tail  for  which  the  ambi- 
polar  solutions,  Eqs.  (B-16)  and  (B-17),  are  applicable.  It  was  found  that 
the  solutions  for  A  =  3000  produced  sufficiently  long  ambipolar  tails  for  all 
the  cases  considered  in  the  present  study  as  shown  in  Fig.  4.  In  this  tail 
region,  where  the  solutions  of  the  two  gegions  overlap,  the  matching  is 
accomplished  by  requiring  that  the  concentrations  and  the  current  densities 
of  the  electrons  and  of  the  ions,  and  the  electric  field  intensity  be  con¬ 
tinuous.  Such  requirements  produce  the  following  algebraic  equations  with 
the  aid  of  the  numerical  solutions  of  Table  1  and  Eqs.  (B-l 5)  and  (30): 


3  ,  _ 
•Pn  +  G 


2<JiG3 


-  - 


G1G2G3 

n/2- 


Jl=° 


(B-18) 


1 


mQ  =  — | 

J, 

G  /°1  -'fc’loV 

+  1 

1 

3\  /J 

(B-19) 


where 


(B-20) 
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^r> 

■'  * 


The  temperature  ratios  T'^/T  appearing  in  Eqs.  (B-16)  through  (B-20), 
and  in  Table  1  are  the  average  temperature  ratios  for  the  sheath.  The 
solution  of  Eqs.  (B-18)  and  (B-19)  gives  r|0  and  mQ  with  which  m.,  m g  and 
E  profiles  can  be  obtained  from  Eqs.  (B-15)  and  (30)  for  the  ambipolar 
region.  The  numerical  results  for  a.,  ae>and  H  obtained  for  the  sheath 
region  can  be  transformed  to  rrn,  mg  and  E  by  the  use  of  t)q  and  and  the 

relationships  between  a  and  m.and  H  and  E,  respectively.  Typical  electron 
and  ion  concentrations  and  electric  field  intensity  profiles  are  shown  in 
Fig.  5.  The  profiles  shown  in  the  figures  are  those  corresponding  to  the 
electron  temperatures  for  Tg^  =  7  20  (Fig.  3). 

The  potential  drop  across  the  boundary  layer  is  obtained  by  integrating  the 
electric  field  intensity  over  the  physical  thickness  of  the  boundary  layer. 

Thus 


Since  E  is  given  in  closed  form  and  pg/p  can  be  obtained  from  Ref.  5,  the 
second  of  these  integrals  may  be  readily  evaluated.  However,  H  is  known 
only  numerically  so  that  the  first  integral  is  more  difficult.  However, 
Chung  has  shown  (Ref.  2)  that 


/p)H  dz 


(B- 22) 
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f) 

where  the  two  integrals  on  the  right-hand  side  have  been  computed  for  sheath 
temperature  ratios  in  the  span  of  interest  (Table  1). 


/ 


/ 
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